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^ . Abstract 

In the final few years of his life, Julian Schwinger proposed that the "dynamical Casimir 
■ effect" might provide the driving force behind the puzzling phenomenon of sonoluminescence. 

. Motivated by that exciting suggestion, I have computed the static Casimir energy of a spher- 

ical cavity in an otherwise uniform material with dielectric constant e and permeability /i. 
As expected the result is divergent; yet a plausible finite answer is extracted, in the leading 
uniform asymptotic approximation. That result gives far too small an energy to account 
for the large burst of photons seen in sonoluminescence. If the divergent result is retained 
(which is different from that guessed by Schwinger), it is of the wrong sign to drive the 
effect. Dispersion does not resolve this contradiction. However, dynamical effects are not yet 
included. 
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^ : 1 Introduction 



Oh! 



In a series of papers in the last three years of his life, Julian Schwinger proposed ^ that the 



D . dynamical Casimir effect could provide the energy that drives the copious production of photons in 
^ ' the puzzling phenomenon of sonoluminescence P, |^. In fact, however, he guessed an approximate 
(static) formula for the Casimir energy of a spherical bubble in water, based on a general, but 
^ ■ incomplete, analysis 0|. He apparently was unaware that I had, at the time I left UCLA, completed 
■ the analysis of the Casimir force for a dielectric ball 0. It is my purpose here to carry out the very 
straightforward calculation for the complementary situation, for a cavity in an infinite dielectric 
medium. In fact, I will consider the general case of spherical region, of radius a, having permittivity 
e' and permeability n', surrounded by an infinite medium of permittivity e and permeability /i. 

Of course, this calculation is not directly relevant to sonoluminescence, which is anything but 
static. It is offered as only a first step, but it should give an idea of the orders of magnitude 
of the energies involved. It is an improvement over the crude estimation used in Attempts 
at dynamical calculations exist [^, but they are subject to possibly serious methodological 
objections. Sonoluminescence aside, this calculation is of interest for its own sake, as one of a 
relatively few nontrivial Casimir calculations with nonplanar boundaries p, ^ |10|, |TT], |T2|, |13|, |1^ 



|T5] . It represents a significant generalization on the calculation of Brevik and Kolbenstvedt 
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who consider the same geometry with fie = fi'e' = 1, a special case, possibly relevant to hadronic 
physics, in which the result is unambiguously finite. 

In the next section we review the Green's dyadic formalism we shall employ, and compute the 
Green's functions in this case for the TE and TM modes. Then, in Section 3, we compute the 
force on the shell from the discontinuity of the stress tensor. The energy is computed similarly in 
Section 4, and the expected relation between stress and energy is found. Estimates in Section 5 
show that the result so constructed, even with physically required subtractions, and including both 
interior and exterior contributions, is divergent, but that if one supplies a plausible contact term, a 
finite result (at least in leading approximation) follows. (Physically, we expect that the divergence 
is regulated by including dispersion.) Numerical estimates of both the divergent and finite terms 
are given in the conclusion, and comparison is made with the calculations of Schwinger. 

2 Green's Dyadic Formulation 



I follow closely the formulation given in |TD|, We start with Maxwell's equations in rationalized 
units, with a polarization source P: (in the following we set c = /i = 1) 

-V X E = , V ■ B = 0, (1) 

dt ' ^ ^ 

where, for an homogeneous, isotropic, nondispersive medium 

D = eE, B = /iH. (2) 

We define a Green's dyadic F by 

E(r, t)= J {dr') dt' r(r, t; r', t') ■ P(r', t') (3) 
and introduce a Fourier transform in time 

du 



r r 



t;r',0= rf^e-(*-*')r(r,r';.;), (4) 



where in the following the uj argument will be suppressed. Maxwell's equations then become 
(which define $) 

V X r = icj$ , V ■ $ = 0, 

Iv X $ = -iujtV , V ■ r' = 0, (5) 
/i 

in which T' = F + 1/e, where 1 includes a spatial delta function. The two solenoidal Green's 
dyadics given here satisfy the following second-order equations: 

(V' + uj^efi)r' = --V X (V X 1), (6) 

e 

(V^ + cj^e/i)$ = iujfiV X 1. (7) 
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They can be expanded in vector spherical harmonics [O, M defined by 
as follows: 

r'(r,r') = {fi{r,v')l^Un) + x (^^(r, r')X,„(^]) 



Irn 



Im 



$(r,r') = ^ ^K^,r')X,„(^]) - - V x /,(r, r')X,„(^]) . 



When these are substituted in Maxwell's equations (|]) we obtain, first, 

91 = 91, fi = fi + -\s{r - r')Ximm, 
and then the second-order equations 

{Di + u^fie)giir, r') = i^/i J dn" X^Xn") ■ V" x 1, 

(A + ^V)/K^,r') = --J dQ"XU^") . V" X (V" X 1), 



e 



where 



These equations can be solved in terms of Green's functions satisfying 

(A + ^'e/i)i^/(^ r') = -^6{r - r'), 

which have the form 

Fi{r,r') ■■ 

where 



ik'ji{k'r<)[hi{k'ry) - Aj,(A;'r>)], r,r' < a, 
ikhi{kr^)[ji{kr<) — Bhi{kr<)], r,r' > a, 



k = \uj\y/JIe, k' = \uj\yfi'e', 
and hi = hf^^ is the spherical Hankel function of the first kind. Specifically, we 

gi{r,r') = -zufiV x Gi{r,r')XU^'), 



where Fi and Gi are Green's functions of the form (IIB) with the constants A and B determined 
by the boundary conditions given below. Given Fi, Gi, the fundamental Green's dyadic is given 

by 



Im 



r'(r,r') = Y.\uj'fiFiir,r')XiUmLm 



-iv X Gi{r,r')XU^)X*^{Q') x V' 



(20) 



Now we consider a sphere of radius a centered at the origin, with properties e', fi' in the interior 
and e, /i outside. Because of the boundary conditions that 



E_L, eEr, Br, 



B 



/i 



(21) 



be continuous at r = a, we find for the constants A and B in the two Green's functions in (^ 



Af 
Bf 
Ac 
Bg 



/q7ei{x')e'i{x) — y/VJjiei{x)e'i{x') 
/e/7s/(x')s^(x) — -/FJlsi{x)s\{x') 



A, 



e' nei{x')e'i{x) 



^ ei{x)e'i{x') 



^FJIsi{x')s'i{x) — y/eJ7si{x)s'i{x') 



(22) 
(23) 
(24) 
(25) 



Here we have introduced x = ka, x' = k'a, the Riccati-Bessel functions 



ei{x) = xhi{x), si{x) = xji{x), 



(26) 



and the denominators 



efi'si[x )ei[x) 



e' fiSi{x')e'i{x) 



e'iJ,Si[x )ei[x), 
efi's'i{x')ei{x), 



(27) 



and have denoted differentiation with respect to the argument by a prime. 

3 Stress on the sphere 

We can calculate the stress (force per unit area) on the sphere by computing the discontinuity of 
the (radial-radial component) of the stress tensor: 



T„{a+) 



:28) 



where 



Trr = \{HEl-El)+^^{Hl-Hl)]). 



(29) 



The vacuum expectation values of the product of field strengths are given directly by the Green's 
dyadics computed in Section 2: 



z(E(r)E(r')) = r(r,r'), 

^(B(r)B(r')) = -^V x r(r,r') x V', 



(30) 
(31) 



where here and in the following we ignore 5 functions because we are interested in the limit as 
r' — i> r. It is then rather immediate to find for the stress on the sphere (the limit t' — >■ t is assumed) 



1 /•°° doj 



-iuj(t-t') 



2ia^ J-oo 2tx 



E 

1=1 



2/ + 1 

A-n 



x< e - e 



-a'F,ia+, a+) + (^^-^ + -g^r—r'j G^{r, r') 

^a^G,{a+,a+) + (^-^^ + -^r^A F,{ry) 
jj \ jj, jjior or I 



r=r'=a+ 



r=r'=a+ 



dy 



21 + 1 d 



2a^ J-oo 2tt At: dx 



where y = ua, 6 = (t — t')/a, and 



(32) 
(33) 



In A; A; = In 



{si{x')ei{x) - s'i{x')ei{x)y - ^\si{x')ei{x) + s'i{x')ei{x)y + constant. (34) 



Here the parameter ^ is 



_ 1 



^4 + 1 



(35) 



This is not yet the answer. We must remove the term which would be present if either medium 
filled all space (the same was done in the case of parallel dielectrics [|19]). The corresponding 
Green's function is 



The resulting stress is 



duj 



,(0) ^ I ik'ji{k'r^)hi{k'ry), r,r' < a 
' \ikji{kr^)hi{kry), r,r' > a 



21 4- 1 



(36) 



a-^ J-oo 2tc 



1=1 



All 



{x'[s'i{x')e'i{x') - ei{x')s'l{x')] - x[s'i{x)ei{x) - ei{x)s'l{x)]} 



The final formula for the stress is obtained by subtracting (37) from (^3|): 



(37) 



2a^ J-oo 2ti Att ^ dx 



1=1 



+2x'[s'i{x')ei{x') — ei{x')si {x')] — 2x[s[{x)ei{x) — ei{x)si (x)]}, 



(St 
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where we have now performed a Euchdean rotation, 

y ^ iy, X ^ ix, t = t — t' ^ 2(0:4 — X4) [5 = (X4 — x'^^/a], 

[tTx 2 / TTX 

si{x) si{x) = ^—Ii+i/2{x), ei{x) ei{x) = —Ki+i/2{x). (39) 

4 Total energy 

In a similar way we can directly calculate the Casimir energy of the configuration, starting from 
the energy density 

U=^Jl±tI!l, (40) 
In terms of the Green's dyadic, the total energy is 



E = J (dr) U 



— / r'^drdVt 



2i 

1 r°° duj 



eTrr(r,r) ^Tr V x r(r,r)x V 



(41) 



2i 



V(2Z + 1) / r'dr 
-cxD 2tc ^^-^ Jo 



X i^2e[Fi{r,r) + Giir,r)] + l|-r + G,](r, r')|,,=,| , (42) 

where there is no explicit appearance of e or fi. (However, the value of k depends on which medium 
we are in.) As in fl^ we can easily show that the total derivative term integrates to zero. We are 
left with 

1 foo fli , °° coo 

^ = Tr/ 7^e---^(2/ + l) / r^dr2k^[Fi{r,r)+Gi{r,r)]. (43) 
2i J~oo Zn Jo 

However, again we should subtract off that contribution which the formalism would give if either 
medium filled all space. That means we should replace Fi and Gi by 

F r = I -'^^'MGji{k'r)ji{k'r'), r,r' < a , . 

' \-ikBF^Ghiikr)hi{kr'), r,r' > a ^ ^ 

so then (||^) says 

duj 

[21 + 1) / 

1=1 ''-'^ 

(45) 

The radial integrals may be done by using the following indefinite integral for any spherical Bessel 
function ji: 

dxx^jfix) = -[{{xjiyf - jiixjiY - xjiixji)"]. (46) 



E = - V(2/ + 1) / —e-''^- / r^dr k'^Ap + AG)]l{k'r) + / r^dr k^Bp + BG)h^{kr) \ 

J-00 2tX Ljo Ja J 
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But we must remember to add the contribution of the total derivative term in (^) which no longer 
vanishes when the replacement (|4^) is made. The result is precisely that expected from the stress 



where the derivative is the naive one, that is, the cutoff 6 has no effect on the derivative. 



5 Asymptotic analysis and numerical results 

The result for the stress ([38| ) is an immediate generalization of that given in j^, and therefore, 
the asymptotic analysis given there can be applied nearly unchanged. The result for the energy 
is new, and seems not to have been recognized earlier. 

We first remark on the special case y/ejl = \^e'fi'. Then x = x' and the energy reduces to 



1 noo °o fl 



(4J 



where 



fx — H 



(49) 



/i + /i' 

If ^ = 1 we recover the case of a perfectly conducting spherical shell, treated in [|l^ , for which E is 
finite. In fact (^9|) is finite for all ^, and if we use the leading uniform asymptotic approximation 
for the Bessel functions we obtain 

u^^e. (50) 



Further analysis of this special case is given by Brevik and Kolbenstvedt |TH] . 

In general, using the uniform asymptotic behavior, with x = uz, v = Z + 1/2, and, for simplicity 
looking at the large z behavior, we have 



E 



1 



27ra y/ejl 



dze'^'^^^z^ln 
dz 



1 + 



e/i 



16^4 y^'ij, 



(51) 



which exhibits a cubic divergence as 5 — * 0. To be more explicit, let us content ourselves with 
with the case when e — 1, e' — 1 are both small and /x = /x' = 1. Then, the leading u term is 



\2 oo 



E 



> U - 



dz e 



iuz5 , 



d 



ie' - e)2 



64a 



16 



dz(l + z^Y 
e' - ef 



256a 



(52) 



Here, the last arguable step is made plausible by noting that since 5 = r/a the divergent term 
represents a contribution to the surface tension on the bubble, which should be cancelled by a 
suitably chosen counter term (contact term). This argument is given somewhat more weight 
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- 1 = . .o/,,2 - (54) 



by the discussion in |2D|. In essence, justification is provided there for the use of zeta-function 
regularization, which directly gives the finite part here: 

32™ 2 64a V 47 ' ^ ^ 

because Y^'iZo ^'^ = (2~* — 1)C(~'5) vanishes at s = 2. 

Alternatively, one could argue that dispersion should be included crudely modelled 

by 

If this rendered the expression for the stress finite [we consider the stress, not the energy, for it is 
not necessary to consider the dispersive factor d{ue{uj))/duj there], we could drop the cutoff 6 and 
the sign of the force would be positive: (at last, we set e' = 1) 

where zq = uoQa/v. As z/ — oo, zq — > 0, and the integral here approaches tizq/2, and so 

(eo - 1)^ - (60 - 1)^ 3 

if we take as the cutoff of the angular momentum sum Vc ~ ^^o^. The corresponding energy is 
obtained by integrating — 47ra^jF, 

which is of the form of (|S2|) with 1/ 5 ujq/ A. 



6 Conclusions 

So finally, what can we say about sonoluminescence? To calibrate our remarks, let us recall (a 
simplified version of) the argument of Schwinger On the basis of a provocative but incomplete 
analysis he argued that a bubble (e' = 1) in water (e ^ (4/3)^) possessed a positive Casimir 
energy^ 

Ec / ^h^-k 1 F 1 F , 58 

3 J (27r )3 2 \ ^e) 12t: \ ^t)' ^ ^ 

■^Inconsistently, for then zq ~ 1. If zq = 1 in (^5|), however, the same angular momentum cutoff gives 5/12 of 
the value in (|5^). 

^Note, for small e — 1, Schwinger's result goes like (e — 1), indicating that he had not removed the "vacuum" 
contribution corresponding to (|37|) . This is the essential physical reason for the discrepancy between his results 
and mine. 
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where i^' is a wavenumber cutoff. Putting in his estimate, a ~ 4 x 10^'^cm, K 2 x lO^cm^^ 
(in the UV), we find a large Casimir energy, Ec ~ 13MeV, and something hke 3 milhon photons 
would be liberated if the bubble collapsed. 

What does our full (albeit static) calculation say? If we beheve the subtracted result, the last 
form in (|52|) , and say that the bubble collapses from an initial radius = 4 x 10~^cm to a final 
radius = 4 x 10~^cm, we find that the change in the Casimir energy is AU ~ +10~^eV. This 
is far to small to account for the observed emission. 

On the other hand, perhaps we should retain the divergent result, and put in reasonable cutoffs. 
If we do so, we have 

E = ^ ' a^K^ ~ -4 X lO^eV, (59) 

perhaps of acceptable magnitude, but of the wrong sign. The same conclusion follows if one uses 
dispersion, as (|57D shows. 

So we are unable to see how the Casimir effect could possibly supply energy relevant to the 
copious emission of light seen in sonoluminescence. Of course, dynamical effects could change 
this conclusion, but elementary arguments suggest that this is impossible unless ultrarelativistic 
velocities are achieved. (See also 0.) Yet the subject of vacuum energy is sufficiently subtle that 
surprises could be in store. A more complete analysis will be provided elsewhere. 
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